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Abstract 

We study non-standard interactions (NSIs) at reactor neutrino experiments, and in particular, 
the mimicking effects on $13. We present generic formulas for oscillation probabilities including 
NSIs from sources and detectors. Instructive mappings between the fundamental leptonic mixing 
parameters and the effective leptonic mixing parameters are established. In addition, NSI correc- 
tions to the mixing angles #13 and 612 are discussed in detailed. Finally, we show that, even for 
a vanishing #13, an oscillation phenomenon may still be observed in future short baseline reactor 
neutrino experiments, such as Double Chooz and Daya Bay, due to the existences of NSIs. 
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I. INTRODUCTION 



Neutrino oscillations have successfully turned into the most plausible description of neu- 
trino flavor transitions. At the moment, the most important task in neutrino physics is to 
accurately determine the neutrino parameters, especially the neutrino mass hierarchy and 
the leptonic mixing angle #13. In this work, we will concentrate on the leptonic mixing 
parameters, and in particular, the parameter # 13 . Recently, so-called non-standard interac- 
tions (NSIs) have been extensively studied in the literature. Such NSIs could affect neutrino 
oscillations in a similar way as normal matter affects them. Therefore, if present, NSIs will 
affect the determination of the fundamental neutrino parameters. 

In this work, we will mainly investigate measurements of the fundamental leptonic mixing 
angles #13 and 612 at reactor neutrino experiments. Since reactor neutrino experiments 
such as the future Double Chooz [l] and Daya Bay [2] experiments as well as the existing 
KamLAND experiment |3j] have relatively short baseline lengths, normal matter effects are 
negligible. This also holds for NSI effects during propagation of neutrinos. Thus, we will 
only assume that the NSI effects are present at sources and detectors. 

In Ref. [l], the basic formalism and different neutrino states for source and detector 
including NSIs (or "new physics") were first presented. Later, NSIs in sources and detectors 
have been discussed using amplitudes that describe the neutrino sources and detectors. Such 
works have been carried out for long baseline neutrino oscillation experiments in general \^ 
as well as for neutrino factories in particular 6]. Recently, a study on the impact of NSIs 
at reactor and accelerator neutrino experiments has been performed jj]. Especially, the 
authors of this work derive first-order series expansions for oscillation probabilities including 
NSIs from sources and detectors. Explicit upper bounds on parameters describing NSIs 
from sources and detectors exist. However, these bounds are only generic and given by 
e s a p = 0(0.1) for NSIs at sources from universality in lepton decays and elg = 0(0.2) for 
NSIs at detectors from universality in pion decays (si, ^] . 

This paper is organized as follows. In Sec. [Ill we will present general formulas for pa- 
rameter mappings between the fundamental leptonic mixing parameters and the effective 
leptonic mixing parameters due to the effects of NSIs, and we will give expressions for os- 
cillation probabilities. Then, in Sec. IHfl we will discuss reactor experiments and how these 
could be influenced by NSIs and what the outcome would be for the mixing angles #13 and 
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#12- Finally, in Sec. IIVI we will summarize our results and present our conclusions. 



II. ANALYTIC FORMALISM 



For a realistic neutrino oscillation experiment, in the presence of non-standard neutrino 
interactions, the neutrino states produced in the source and observed at the detector can be 
treated as superpositions of pure orthonormal flavor states: 



V l3=e,fi,T 
P \ a=e,ii,T 



(1) 



(2) 



where the superscripts 's' and 'of denote the source and the detector, respectively, with the 
normalization factors being given by 1 



[(l + e s ) (l + £ s t)] 



act ' 



(3) 
(4) 



Note that the states and (u d \ are no longer orthonormal states because of NSIs. Since 
different physical processes take place at the source the and detector, the NSI parameter 
matrices e s and e d are arbitrary and non-unitary in general. In the minimal unitarity 
violation model (MUV) 



matrix 



16 



3, 



ii 



12 



13l . Il4 , 1 la ], where the unitarity of the leptonic mixing 



171 ] is slightly violated by possible new physics effects, the non-unitary effect can 



be regarded as one type of NSIs with the requirement e s = e d ^ ? 

Since in a terrestrial neutrino oscillation experiment, the Earth matter effects [18|, H9( are 
more or less involved, the propagation of neutrino flavor states in matter is governed by the 
effective Hamiltonian 



H = H + H m + if NSI 

= — ZJdiagK, ml ml)U^ + diag(V cc , 0, 0) + V cc e r 



(5) 



1 Note that, in calculating the number of events, the normalization factors are canceled with the NSI factors 
in charged-current cross-sections. However, for running short baseline reactor neutrino experiments, the 
neutrino fluxes are directly measured by using a near detector, and not a Monte Carlo simulation. Hence, 
the normalization factors should be taken into account. (See also Ref. 10( for a detailed discussion.) 

2 In the MUV model, a neutral-current contribution cannot, in principle, be rewritten as a global phase in 
the oscillation amplitude, and thus, it affects the oscillation process. 
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where Vcc = V^GpN e arises from coherent forward scattering and N e denotes the electron 
number density along the neutrino trajectory in the Earth. Different from e s and e d , e m is 
an exact Hermitian matrix describing NSIs in matter, and its current experimental bounds 



can be found in Refs. 



21 



221 ] . Here the superscript 'm' is used in order to distinguish 



NSI effects in the Earth matter from those in neutrino sources and detectors. The vacuum 
leptonic mixing matrix U is usually par ametrized in the standard form by using three mixing 



angles and one CP violating phase 



23| 



u = o 23 u s o 13 u}o 
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S 12 S 23 S 13 e 



s 13 e 



-A 
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-C-inS 
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s 12 c 23 s 13 e CqoC 



23^13 



(6) 



/ 



where Us = diag(l, 1, e lS ), and Oij is the orthogonal rotation matrix in the i,j plane with 
c i:( - = cos 9^ and = sin 6^ (for ij = 12, 13 and 23). In analogy to the vacuum Hamiltonian 
Hq in Eq. (jSJ), the effective Hamiltonian in matter H can also be diagonalized through a 
unitary transformation 



// ^rdiag (rh\ ~ 2 2 



, m 2 , m 3 



(7) 



where m? 



1,2,3) denote the effective mass squared eigenvalues of neutrinos and U is 



the effective leptonic mixing matrix in matter. 

Now, we include all the NSI effects into the oscillation processes, and arrive at the am- 



plitude for the process v s a 
1 



iHL I , ,s \ 



(1 + e d ) pP A lp (1 + e s ) 



N' a N* 



(l + eyA 1 (! + £'' 



j3a 



N' a N* 



[A + e s A + Ae d + e s Ae d ] , (8) 



where L is the propagation distance and the explicit form of A is a coherent sum over the 
contributions of all the mass eigenstates 



A 



E.m-L 



(9) 
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Inserting Eq. (jHJ) into Eq. one can directly obtain 

1 



A aP (L) 



N»N* 



7,i 



w ■» m ■ L ^ — , >, a T 

+ E ^^e^e-'^r + £ e^U^U^- 



7,1 



7>P>« 



— F 



N' a N$ 



0i 



+ E/ ^ifiafi^i 



e cey £ pfiU ^iU pi 



e 2E 



7 7,P 

In order to compare Eq. (fTUj) with the standard oscillation amplitude given in Eq. 
rewrite A^L) as 



with 



a/3 



(10) 



), we 



(11) 



(12) 



It can be clearly seen that only the ath row of e s and the /3th column of e d are relevant to 
the transition amplitude. In the limit e — > 0, Eq. fTTTj) is reduced to the standard oscillation 
amplitude in matter. 

With the definitions above, the oscillation probability is given by 



4) = \A a p{L)\ 2 

= E - 4 E mw;) sm 2 Al%lL 

-2E Im (^7$)sin 

i>3 



AE 



2£ 



(13) 



A salient feature of Eq. fflBl is that, when a ^ (3, the first term in Eq. ffl3|) is, in general, not 
vanishing, and therefore, a flavor transition would already happen at the source even before 
the oscillation process and is known as the zero-distance effect 24J. Although the effective 
mixing matrix in matter U is still unitary, the presences of NSIs in the source and detector 
prevent us from defining a unique CP invariant quantity like the standard Jarlskog invariant 



251 ] . New CP non-conservation terms, which are proportional to the NSI parameters and 



have different dependences on L/E, will appear in the oscillation probability. Another 
peculiar feature in the survival probability is that, in the case of a = (3, CP violating terms 
in the last line of Eq. (I13p should, in principle, not vanish. Note that Eq. (I13p is also valid 
in the MUV model and could be very instructive for analyzing the CP violating effects in 
the MUV model in future long baseline experiments. 

III. REACTOR NEUTRINO EXPERIMENTS 

Reactor neutrino experiments with short or medium baselines are only sensitive to the 
survival probability P{p s e — > vf). The typical energy of antineutrinos produced in nuclear 
reactors is around a few MeV, which indicates that the Earth matter effects are extremely 
small and can safely be neglected. Hence, we take (U ~ U , rhi ~ m«) or effectively set 
Vcc — in Eq. (jSJ). As mentioned above, among all the NSI parameters, only e s ea and e^ e 
are relevant to our discussion. It has been pointed out that for realistic reactor neutrino 
experiments, the leading-order NSIs are of the V ± A type, and the relation e s ea = e d * e holds 
well [7]. Therefore, we assume e s ea = e d * e = \ee a \^ ea in the current consideration and neglect 
the superscript V throughout the following parts of this work. It can be seen from Eq. f[T2"j) 
that the imaginary parts of the parameters Jl e disappear, and hence, the corresponding 
i>l — > z/g oscillation is a CP conserved process. 

Similar to the case without NSIs, one may define the effective mixing angles #13 and 
612, in which all the NSI effects are included. For the smallest mixing angle #13, we take 
a = j3 = e and i = 3 in Eq. (TT21) together with the standard parametrization defined by 
Eq. (jSJ), and obtain the mapping between #13 and #13 

S 13 = S 13 + 2Sl 3 Cl3 [S23 OOS(5 - e/ J|e e/J | + C 23 C0S(5 - (per)\£er\ 

S 2 3 C0S(5 - (pee - 4>e^)\See\\£e^\ ~ C 2 3 COs(5 - ee - er ) \e ee \ \e eT |] 

, / 2 2 2 \\ |2 , / 2 2 2 \ I |2 

+ l s 23 C 13 — s 13)\ £ e^\ + l c 23 C 13 — S l3)\ £ er\ 

+2s 2 3C 2 3Ci3 COS(0 

efj, 0er) |^e/^| |^er I + 0(e 3 ) , (14) 

where the third-order terms in e are neglected. As for the effective mixing angle (9 12 , we take 
a — (3 — e and i = 2, and obtain 

■812^3 = s i2 c i3 + 2si2Ci 2 Ci 3 [c 23 cos(0 eM ) \e e/l \ - s 23 cos(0 er ) |e eT |] 

-2s 2 12 s 13 c 13 [s 23 cos(5 - eM )|e eA1 | + c 23 cos(5 - eT )|e eT |] + 0{e 2 ) . (15) 
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Since the NSI parameters should not be comparable to the sizable mixing angle #12, only 
the first-order terms in e are taken into account in Eq. (j!5p . 

With the help of the effective mixing angles #13 and #12, the survival probability reads 



P(pl^pj) = 1 - cos 4 13 sin 2 20 12 sin 2 - cos 2 § 12 sin 2 20 13 sin 2 ^ L 

AE AE 



sin 2 12 sin 2 26 13 sin 2 - . (16) 



A. Short baseline reactor experiments and #13 

The forthcoming two improved short baseline reactor neutrino experiments Double Chooz 
and Daya Bay are planned with the same goal of searching for the smallest leptonic mixing 
angle 13 . Both of these two experiments make use of the same concept: one near detector 
is placed a few hundred meters from the core of the nuclear power plant in order to reduce 
systematic errors and one far detector is located at distance (L ~ 1 — 2 km) close to the 
first maximum of the survival probability caused by the large mass squared difference Am 2 x . 
The v e — > v e channel is dominated by the atmospheric oscillation dip, which allows us to 
safely neglect the term containing Am^ in Eq. (flBT) . and we arrive at 

Piyl - vt) * 1 - sin 2 2~9 l3 sin 2 . (17) 

Since NSIs are only sub-leading order effects, higher-order terms proportional to es1 3 can 
be ignored, and then the effective mixing angle #13 in Eq. (THI) approximates to 

5 i3 = s i3 + 2s i3 [S23 cos(5 - 6(U ) |e e/i | + c 23 cos(5 - er ) |e eT |] 

) + 0(e\esl 3 ) . (18) 

Note that s 2 3 is invariant with respect to the exchange e efl <-> e eT , and obviously, in the limit 
e — > 0, #i3 equals 6*i 3 . Equation (flBl) clearly shows how the mixing angle 9 13 is modified by 
NSIs. Some comments are in order: 

• The contributions coming from the NSI parameter e ee are always correlated with 
higher-order corrections, and hence cannot be well constrained in a reactor experiment. 
However, it induces an enhancement of the total neutrino flux at a near detector, 
which appears as an overall factor in the oscillation probability if we do not normalize 
neutrino states as in Eq. ([IT) . Due to the flux uncertainty in reactor experiments, it is 



very hard for this enhancement to be observed [7] 



For a given set of NSI parameters, sin #13 is a quadratic function of sin #13. Thus, there 
exists a minimum of sin 2 # 13 at the position 



Sl 3 |min = -S 23 C0S(5 - e u)keu| ~ C 23 COs(5 - <f) eT )\e 



(19) 



and the minimum value of sin 2 # 13 is given by 



s 2 I • 

°13 |min 



^23l^e/i| SID. 2 



23l c er| SID 2 



+S23C23| I sin(<5 - 6(U ) sin(<5 



^23 I ^e/x I "I - C 23 l^e-r | 
+2|^ e/ i| |£ eT |s23 c 23 COS(0 

e/i 0er) 



for Si 3 | min > 



• (20) 



for Sl 3 |min < 



Since the fundamental 6*13 cannot be well distinguished from the effective #13 measured 
in an experiment, the mimicking effects of NSIs play a very important role in the 
small 6*13 region. Even if the true value of 9 13 is too tiny to be detected, we may 
still hope to obtain an oscillation phenomenon in reactor experiments. On the other 
hand, compared to 6*13, #13 may also be remarkably suppressed by NSIs, which makes 
the current experiments quite pessimistic. Note that mimicking (or "fake") values 
of 6>i3 due to so-called damping effects have been investigated in Ref. |2a|. Such 
damping effects could arise from decoherence-like damping signatures (e.g. wave-packet 
decoherence related to production and detection processes). Thus, damping could fake 
values of #13, and therefore the value of #13 would turn out to be smaller than one 
expects. 

We illustrate the mappings between #13, e, and 6*13 in Fig. [TJ In our numerical cal- 
culations, we use the exact analytical formulas and do not make any approximations. 
We also adopt the central values of other relevant parameters from the global fit given 



in Ref. [23]. Without loss of generality, we take \e 



e/x I 



|e| in our analysis , 



and allow all the CP violating phases to vary from to 2ir. For a given value of #13, 
which is in fact the parameter measured in experiments, the true values of #13 may 
be remarkably different, i.e., there exists a degeneracy in #13. Therefore, one has to 



Since we do not make any constraint on the CP violating phases, the numerical results would almost be 



the same for the case 
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N Ml 

FIG. 1: Mappings between #13, #13, and NSI parameters e a p. In the left plot, we assume < 
#13 < 10° as experimental constraints. The shaded areas correspond to different upper bounds on 
the effective mixing angles. For the right-hand plot, \e\ is allowed to vary from to 0.05, and the 
gray shadings represent |e| < 0.05, \e\ < 0.01, and |e| < 0.001, respectively, with darker regions for 
smaller |e|. All CP violating phases are treated as free parameters and allowed to vary from to 
2vr. 



disentangle the parameter #13 from the NSI parameters. Within the present upper 
bound #13 < 10° [281 ] . #13 may approach 14° at large e regions. In the case #13 < 5°, 
there is still a widely allowed range 1° < #13 < 8° with respect to a large e. Even 
if #13 is too small to be measured in a reactor experiment, i.e., #13 < 3°, a discovery 
search of a non-vanishing #13 may still be carried out at future neutrino factories, 
where the source of neutrinos is a muon storage ring with very clean muon decay and 



quite limited room for NSIs 



29]. 



• In Fig. [21 we show the oscillation probabilities with respect to the NSI parameters. 
The upper plot in Fig. [2] indicates that mimicking oscillation effects, which are induced 
by sizable NSIs, can be observed in despite of a negligible #13. Once other type of 
neutrino oscillation experiments can help us to fix the true value of #13, the mimicking 
effects will provide us with the opportunity to search for NSIs in neutrino sources and 
detectors. 

The oscillation process expressed in Eq. (fTTl) is actually CP conserved. However, the CP 
violating phase 5 enters the oscillation probability explicitly, and so does the leptonic mixing 
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FIG. 2: NSI corrections to the oscillation probabilities P{p% — > vf) in a short baseline experiment. 
The shadings correspond e < 0.05, e < 0.01, and e < 0.001, respectively. The values of #13 are 
labeled on the plots. For other mixing parameters, we use the central values given in Ref. j^ . 
Here we take the average energy of reactor neutrinos E = 3 MeV. 
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angle 823. It is then very helpful to extract information on leptonic CP violation and 823 
by analyzing the corresponding disappearance channel together with future long-baseline 
appearance experiments. 



B. Medium baseline reactor experiments and #12 

The current medium baseline reactor neutrino experiment KamLAND receives v e from 
nuclear reactors located at an average distance L ~ 180 km. In order to improve the accuracy 
of current measurements, the next generation experiments should take the baseline length 
of about 50 km, which is close to the first minimum of the survival probability related with 
the small mass squared difference Am^. 

In neglecting contributions from # 13 , the corresponding oscillation probability reads 

P(pt - t) - 1 - sin 2 26 l2 sin 2 , (21) 

where the effective mixing angle 612 approximates to 

*12 = s 12 + 2Si 2 Ci2 [c 23 COs(0 eM )|e eM | — S 2 3 COs(</> 

er) l^er 

\}+0(es 13 ,s 2 13 ) . (22) 
We now discuss how the NSIs affect the leptonic mixing angle #i 2 : 

• Similarly, there is no contribution coming from e ee at leading order. Thus, we can 
observe that reactor experiments are not sensitive to e ee . 

• Compared to the tiny #13, the magnitude of 8u is more sizable. Hence the NSI effects 
cannot mimic an effective 812 with a vanishing 612- However, NSI effects may dramat- 
ically modify the observed mixing angle 812. We plot # 12 as a function of 812 and e in 
Fig. [3j In the large e re gion s, the true value of 812 may be close to the bi- maximal 



mixing value 45° 



30 



31 



32|. On the other hand, the lower bound # 12 > 26° deviates 



much from its tri-bimaximal mixing pattern 33j, |34j]. Figure |3] indicates that, in the 
presence of NSIs, even if 8yi can be well measured, there is still a large room of 812 for 
various flavor symmetric models. 

The oscillation probabilities of medium baseline reactor experiments are illustrated in 
Fig. HI We take the best-fit values of 813 and 823 in our numerical calculations j^ . 
Hence, the oscillation behavior around L ~ is mainly induced by Am^. It can be 
clearly seen that NSI corrections are more significant for a smaller # 12 . 
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FIG. 3: Mappings between #12, #12, and NSI parameters e a p . The shaded areas in the left-hand 
plot correspond to different upper bounds on 6\2- For the right-hand plot, |e| is allowed to vary 
from to 0.05, and the gray shadings represent |e| < 0.05, |e| < 0.01, and |e| < 0.001, respectively, 
with darker regions for smaller |e|. As in Fig. [TJ we allow all the CP violating phases to vary from 
to 2vr. 



Unlike the mapping between # 13 and #13, only 823 is entangled in Eq. fl22l) . Hence, we 
may also acquire useful constraints on 623 through precision measurements of 6*12 and NSI 
parameters in future experiments. 



C. Correlations between #13 and O12 

A crucial question for future experiments is how to distinguish real mixing parameters 
from NSI effects. As discussed above, there is no hope to extract the fundamental mixing 
angles from a single reactor neutrino experiment. Since the NSI effects bring in intrinsic cor- 
relations between the effective mixing parameters, a combined analysis of both appearance 
and disappearance channels should be able to help us to determine NSI effects. 

In Fig. [SJ we show the confidence regions of #13 and #12, constrained by the current global 
fit of neutrino oscillation data j^]. Because of the experimental uncertainties associated 
with #i3 and #12, the allowed parameter spaces for #13 and 612 are quite wide. The true value 



of 6x3 can achieve the range of the Cabibbo angle within 



some light on the quark-lepton complementary models 



32 



a confidence level, which shades 
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361 ] . We want to stress that 
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FIG. 4: NSI corrections to the oscillation probabilities P{y s e — > v~) in a medium baseline experi- 
ments. The shadings correspond e < 0.05, e < 0.01, and e < 0.001, respectively. The values of #12 
are labeled on the plots. The other input parameters are the same as in Fig. [2j 
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FIG. 5: The la, 2a, and 3a confidence regions of fundamental neutrino mixing angles #13 and #12, 
constrained by the current global fit of neutrino oscillation data. The value |e| = 0.05 is assumed 
with all the CP violating phases being allowed to vary from to 2ir. 

our computations depend on the input parameter e and can just serve as a rough illustration. 
D. A discussion on #23 and 8 

Finally, we briefly discuss the NSI corrections to #23 and 5. One may define an effective 
mixing angle #23 by using the analogous way that we performed above for # 13 and (9 12 . How- 
ever, since reactor neutrino experiments are only sensitive to the first row of the leptonic 
mixing matrix, the effective 623 loses its meaning. In future long baseline /3-beam experi- 
ments or neutrino factories, where different types of NSIs are involved in the production, 
propagation, and detection processes, one cannot simply employ the language of effective 
mixing parameters as in reactor neutrino experiments. However, the generic formulas given 
in Eq. (113j) are still valid and very helpful for us in order to figure out NSI effects. A de- 
tailed and joint numerical analysis based on Eq. (fT3|) should be very meaningful and will be 
elaborated elsewhere. 
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IV. SUMMARY 



In this work, we have studied NSI effects in reactor neutrino experiments, and in particu- 
lar, the mimicking effects on 13 . We first presented the most general formulas of oscillation 
probabilities with all NSI effects at production, propagation, and detection processes being 
considered. Instead of directly discussing oscillation probabilities, we took use of a more 
straightforward method, which started from the effective amplitude and derived instructive 
mappings between fundamental mixing angles (6*13, #12) and effective NSI corrected mixing 
angles (#13, 9i 2 ) in reactor neutrino experiments. The analytical relations clearly show how 
these mixing angles are affected by NSIs. We have also illustrated the NSI effects at short 
and medium baseline reactor experiments. We found that the mixing angles measured in 
reactor neutrino experiments could be dramatically modified by NSIs at the neutrino source 
and detector. The mimicking effects induced by NSIs play a very important role in a short 
baseline experiment, especially in the case of a tiny 13 . Even for a vanishing 13 , the forth- 
coming Double Chooz and Daya Bay experiments could still perform a discovery search of 
an oscillation phenomenon, which should totally be attributed to NSI effects. 

From the phenomenological point of view, two different and complementary oscillation 
experiments are needed in order to constrain corresponding NSIs. The measurement of NSI 
parameters should be one of the most interesting topics of experimental physics in future 
even before the discovery of leptonic CP violation. 
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